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Abstract 

We study factorizations of rational matrix functions with simple poles on the Riemann sphere. For 
the quadratic case (two poles) we show, using multiplicative representations of such matrix functions, 
that a good coordinate system on this space is given by a mix of residue eigenvectors of the matrix and its 
inverse. Our approach is motivated by the theory of discrete isomonodromic transformations and their 
relationship with difference Painleve equations. In particular, in these coordinates, basic isomonodromic 
transformations take the form of the discrete Euler-Lagrange equations. Secondly we show that dPV 
equations, previously obtained in this context by D. Arinkin and A. Borodin, can be understood as simple 
relationships between the residues of such matrices and their inverses. 

1 Introduction 

In this paper we are concerned with various ways of introducing coordinates on the space of rational matrix 
functions L(2;) (with simple poles) on the Riemann sphere. Examples like this play an important role in 
various applications, like Yang-Baxter maps and matrix solitons [Ves03, Ves07], Lax equations and isomon- 
odromy transformations on algebraic curves, [Kri02a, Kri02b], and others. Our main interest is related to 
the study of discrete isomonodromic transformations and their relationship with discrete Painleve equations. 
The theory of discrete isomonodromic transformations was developed by A. Borodin for polynomial matrices 
in [Bor04] and adapted to the notion of local monodromy using rational matrix functions by I. Krichever 
in [Kri04]. In [AB06] D. Arinkin and A. Borodin used the theory of d-connections on vector bundles to 
explained the appearance of difference Painleve equations, considered from the geonu^tric point of view of 
H. Sakai [SakOl], in the theory of discrete isomonodromic transformations and later, in [AB07], they intro- 
duced the notion of the r-function of such transformation. These r-functions appear as the gap probabilities 
in the discrete probabilistic models of random matrix type. 

I. Krichever conjectured that discrete isomonodromic transformations can be written in the Lagrangian 
form and that they should be related to the universal symplectic form of Krichever-Phong, [KP98]. In 
[Dzh08], using the methods of [MV91], we verified this conjecture for the quadratic (two-pole) case by using 
the multiplicative coordinates on the space of these matrix functions and finding the explicit formula for the 
Lagrangian. Recently F. Soloviev showed that our Lagrangian symplectic form coincides with the reduction 
of the quadratic symplectic form of Krichever-Phong to some symplectic leaves [Sol09] . 

Unfortunately, the multiplicative coordinates are not easily obtained from other characteristic properties 
of L(2:), such as its residue matrices. This is an obstacle to the generalization of our results to the higher-order 
case. In this paper we argue that one way around this difficulty is to consider, in addition to the residues of 



A. Dzhamay 



Factorization of Matrix Functions 



2 



L(^), the residues of the inverse matrix Jj~^{z). Then good coordinates on the space of L(2;), again in the 
quadratic case, arc given by half of the residues data of L(z) and half of the residue data of L~^(z). Such 
residues also explain the almost symmetric form of the expressions for the multiplicative coordinates of L(z) 
and allow us to recognize dPV equations as simple relations between the residues of L(2;)='=^ and the residues 
of the transformed matrices L(z)^^. 

The paper is organized as follows. In Section 2 we give a short overview of the additive representation of 
L(^;). In Section 3 we study the multiplicative representation of L(2;) and establish the relationship between 
the eigenvectors of the residues of L{z)^^ and the eigenvectors of the left and right-divisors of L(z). We 
also significantly simplify many of the arguments and formulas of [Dzh08] . Finally, in Section 4 we restrict 
our attention to the rank-two case. It this case it is possible to introduce the so-called spectral coordinates 
on the space of L(z) in such a way that the equations relating the spectral coordinates of L(z) and L(z) 
are precisely the difference Painleve equations. Our main observation here is the following. In addition to 
spectral coordinates, divisor of zeroes and poles, and some asymptotic behavior at infinity, the entries of 
the matrix L(z) also depend on a choice of a gauge with respect the action of the group of constant non- 
degenerate diagonal matrices. Understanding the change of this gauge from L{z) to L~^{z) and L(2;) allows 
us to easily obtain the expressions of the multiplicative representation of L(2;) in the spectral coordinates 
and the dPV equations. 



2 Additive form of Rational Matrix Functions 

Let Ij{z) be a rational matrix function on the Riemann sphere, rankL(z) = rn. Let zq be a normalization 
point at which L(z) is regular, lim^_>2o L(z) = Lq, detLo 0, and all eigenvalues of Lq are distinct (we 
assume everything is generic). We can then do a gauge transformation to make Lq diagonal, thus reducing 
the global gauge group to the group of diagonal matrices. Without any loss of generality we can assume 
that zo = oo, and so 

lim L{z) = Lo = diag{/9i, . . .,pm}- (2.1) 

z — *oo 

Next, consider the pole structure of L(z). First, we require that L(z) has only simple poles, denote 
them hy zi, . . . , Zk- Next, we require that the determinant det L(^;), which is a meromorphic function on the 
Riemann sphere that is holomorphic away from zi, . . . , z^, also has only simple poles at those points. This 
condition is equivalent to requiring that = reSx^ L(2:) is a matrix of rank one. Finally, we require that the 
inverse matrix M(2;) = L(z)~^ satisfies similar conditions. To summarize, we have: 

L{z) = Lq -I- ' , where Lo = diag{pi, ... and Lj = ajbj, (2.2) 

detL(z)=pi---p„g|=^4^^, (2.3) 
^ JVl 

L(z)-i = M(z)=Mo-^— ^, where Mo = diag{l/pi,...,l/p„}, M^ = c^dj, (2.4) 

i=i ^ 

and the negative signs in M(z) are for future convenience. We call the above representations of L(z) and 
M.{z) additive representations and the vectors aj, b- (resp. Cj, d-) additive eigenvectors of L(2;) (resp. 
M(z)). By the divisor of L(z) we mean the the divisor T> of its determinant, V = z, — Q. We denote 
the space of matrices L(z) satisfying the conditions (2.2)-(2.4) by . 

An important question is how to choose a good coordinate system on the space M.^. For example, given 
Li, we can determine and bj only up to a common scaling factor. This factor has to be adjusted to ensure 
that detL(z) has zeros at which a complicated condition on tr(Lj). Same problem is present for the 
collection {cj, dj}. Some insight for a good choice of coordinates is provided by the study of the isospectral 
and isomonodromic transformations, which suggests that half of the coordinates should be taken from the 
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residues of L(2) and half from the residues of M(2;). In fact, in the quadratic case, we have the following 
result. 

Theorem 2.1. The vectors (02, d|; a2, b|), considered up to rescaling (i.e., as points inW~^ ) are coordinates 
on the space A4f . To recover 'L'^^{z), consider the function 

'C((x2,x|),(y2,yl)) = (^;2-2:i)log(x|LoX2) + (2i-C2)log(ylx2) + (C2-Ci)log(ylLo V2) + (Ci-^2)log(x|y2). 

(2.5) 

Then 

ai = -^((c2,dl),(a2,bl)); = — ((c2,dl), (a2,bl)); (2.6) 

fir 

ci= _((c2,dl),(a2,bl)); dt = - — ((C2,dl),(a2,bl)). (2.7) 

The proof of this theorem is based on the multiplicative representations of L(2;), which we consider next. 

3 Multiplicative form of Rational Matrix Functions 
3.1 Elementciry Divisors 

To change from the additive to a multiplicative representation, we first need to define the building blocks 
for it. We call such building blocks the elementary divisors. 

Definition 3.1. An elementary divisor is a rational m x m-matrix function B(z) on the Riemann sphere 
of the form 

B(2;) = I H where G = pq^ is a matrix of rank one. (3-1) 

A direct calculation establishes the following elementary facts. 
Lemma 3.1. Let B{z) he an elementary divisor. Then 

(i) detB(2) = (-2 - C,o)/{z - zq), where Co = 2:0 - q^P; 

(ii) B(2)-i = I-G/(^-Co). 

In fact, for us it will be more convenient to fix the points Co and zq on CP^. Thus, we say that a pair 
(Co,-2o) corresponds to an elementary divisor B(z) with the determinant dctB(z) = (z — Q^)l{z — Zq). In 
fact, any generic matrix B(z) with such determinant and normalized by the condition B(z) ^ I as z ^ 00 is 
of the form B(2:) = I + G/\z — zq) with G of rank one and det B(Co) = (more carefully, 'R{z)~^ should be 
regular at zo and B(z) should be regular at Co)- Thus, B(Co) has a left null- vector q''' and a right null- vector 
p. We then immediately get that G = pq^, where we need to normalize the vectors p and q^ so that 
q^p = 2:0 — Co- In a more invariant form this can be written as 

B(z)=I+^5^P?i. (3-2) 
z — Zo qTp 

Prom that point of view, the formula for the inverse matrix follows from the vanishing of the residue of the 
identity B(z)B(2)-i = I = B(z)-iB(z) at zo, 

=i+T — ^-f;:- (3-3) 

The following easy properties of elementary divisors are very useful for what follows- 
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Lemma 3.2. Let 

z- Zo qTp 

Then 
(i) 

BWp=(J^)p <^nd qtB(z)=(^^)qt. (3.4) 

(a) Suppose that at some point z* we have B(2;*)w = v for some vectors v and w. Then 

B(z) = I+^— ((zo-z*)^ + (z*-Co)^) (3.5) 

(i.e., we can determine p from v, w, and ). Similarly, i/w^B(^;*) = v^, 

B(^) = I + — ((zo - ^*)^ + {z* - Co)^) . (3.6) 
2; — ^0 V w P P / 

Proof. Part (i) is immediate. To establish the first formula in part (ii), we solve B(2;*)w = v for p/(q^p) 
and then use (i): 

p zq- z* w z* - Zq w zq- z* 0* - Co V 



qtp Zo - Co q'l'w Zo - Co q^w Zo - Co q^w Zo - Co q^v ■ 
The second formula is proved in the similar way. □ 

Finally, we need the following notation. 
<i Notation: We define a twisting of an elementary divisor B(z) by a constant non-degenerate matrix A to 
be a new elementary divisor ■*-B(z) such that ■*'B(z)A = AB(z), i.e., 



-B(z) = AB(z)A- = I + ^o-Co(Ap)(qtA-^) 

z-zo (qTA-i)(Ap) 



3.2 Factors and Divisors 

We begin with the following important remark. For the additive representation of L(z), the ordering of 
zeroes and poles of dot L(z) is not important, but for any multiplicative representation choosing such an 
ordering is crucial. Thus, from now on our labeling will reflect the fact that (Cs, 2s)-pah corresponds to some 
elementary divisor in L(z). There are two ways to look at the multiplicative structure of 'L{z) — we can 
look at factors or at divisors. 

Definition 3.2. We say that elementary divisors Cs(z) = I + ii-s/{z — Zg) corresponding to pairs (Csj-^s) 
with Hg = mgnj are the factors of L(z) if 

L(z) = LoCi(z)---Cfc(z). (3.7) 

Definition 3.3. We say that elementary divisors B^(z) (resp. B^(z)) corresponding to pairs {Cs^Zs) are 
right (rcsp. left) divisors of L(z) if L(z) = L^(z)BJ(z) (resp. L(z) = Bj^(z)L^(z)) where Lf(z) (resp. L^(z)) 
is regular at Zg. 

The main advantage of the divisors as opposed to the factors is that they can be written explicitly in 
terms of the residues of L(z) and M(z). Note that Cfe(z) = B^(z) and ^"Ci(z) — B'j(z). In particular, in 
the quadratic case fc = 2, there is no essential difference between divisors and factors. 
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Lemma 3.3. Let = res^^ L(z) = a^bj and Mg = — res^^ M(z) = CgdJ. T/ien 

B:(.)=I+^^ (3.8) 
Bi{z)=I+'^^. (3.9) 

Z Zs Us3-s 

Proof. To obtain the formula for the right divisor, we take the residue of L(^;) = LJ(2:)BJ(2;) at Zg to get 
^ (cCV- Then wc take the residue of M(z) = {B'^^izyj-^LKz))-^ at Cs to get c^ ~ pj. The formula 
then follows. The expression for the left divisors is obtained in the same way. □ 

3.3 Refactorization Transformations 

As shown in [Bor04], [Kri04], isomonodromic transformations on the space of rational matrix functions L(^) 
have the form 

L(^) ^ t{z) = R{z + l)L{z)R{z)-\ 
where R(^) is a certain rational matrix function. A simpler version of this axe isospectral transformations 

L(z) ^ L{z) = R(z)L(z)R(z)-i. 

Such transformations arc generated by elementary transformations given by the interc;lianging of factors in 
the multiplicative representation of L(z). Let us now restrict our attention to the quadratic case and consider 
this transformations with Ii{z) = Bj'(2;): 

L(2) = B^(^)LoB^(z) ^ t{z) = B^(z)B^(z)Lo = B^(z)LoB^(z). (3.10) 

Note that the isospectral transformation fixes the divisor 2?, and the isomonodromic transformations shift 
T> by an integer vector. 

To study these transformation from the point of view of the discrete Euler-Lagrange equations it is 
necessary to construct a configuration space Q, a Lagrangian function £ G ^(Q x Q), and a map t] : QxQ^ 
A4^ such that r] maps the Lagrangian dynamics to the refactorization dynamics. That is, denoting by Q the 

previous and by Q the next point of the discrete dynamics, we want J?(Q,Q) = L(^;) and r?(Q,Q) = L(2;), 
see [MV91], [VesQl] for the general description of this approach, and [Dzh08] for our specific situation. 
Writing 

t{z) = B^(^)B^(^)Lo = B[{z)LoB^2{z) (3.11) 

and using the language of divisors and Lemma 3.3, we see that ci = ai, b| = d|, a2 = L0C2, and d^Lo = bj. 

Thus, if we put Q — (a2,b|), Q = (a2,b|) = (LoC2,d|), the statement of Theorem 2.1 follows from 

Theorem 3.1 of [Dzh08]. Below, we restate that theorem in these coordinates (we also improve those formulas 
by using Lo instead of its root) . 

Theorem 3.4. 

(i) The map r] : Q x Q M.^ is given by 

m,Q) = Hz) = fi+ — ({z. - Ci)S^ + (Ci - (2)"^''''^" 

\ z- Z2 y h[a2 

V ^-^1 V b{Lo^a2 




(^1 - C2)rf^ + (C2 - Ci)7T;r^ I I ; (3-12) 
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(ii) The equations of motion (Q,Q) = $(Q,Q) of both the isospectral and isomonodromic dynamics in 
these coordinates are given by the discrete Euler-Lagrange equations with the Lagrangian function 

£(X, Y, t) = {Z2 - zi{t)) log(xlx2) + (01 (i) - C2) log(ylLo-ix2) (3.13) 
+ (C2 - Ci(i)) log(ylLo V2) + (Ci(i) - ^2) log(xly2), (3.14) 

where in the isomonodromic case zi{t) = z\—t, (i{t) = (i — t, and in the isospectral case zi{t) = z\, 
C,i{t) = Ci and C{X,Y) is time-independent. 

The proof is based on computing the residues of L(z) at the points Zi and residues of L(z)~^ at the 
points Ci, followed by the application of Lemma 3.2(ii) to the row and column parts of the resulting rank-one 
matrices. 

4 Isomonodromic Transformations and dPV 

In [AB06], Arinkin and Borodin showed that for rank-two matrices isomonodromic transformations above, 
when written in a special coordinate system, are given by the difference Painleve equations. In this section, 
choosing the dPV case as an example, we show that these equations appear explicitly as relations between 
the residues of L='=^(z) and L='=^(z). 

4.1 Spectral Coordinates 

In the quadratic (two-pole) case the space of the rank-two matrices L(z) satisfying the requirements (2.2-2.3) 
can be described using different parameters. These parameters come in two groups. The first group, that 
we call the type of 'L{z), consists of the zeroes and poles of the determinant of 'L{z) and some asymptotic 
data at z = 00. The space of L(2;) of the fixed type is two-dimensional, and the sec;ond group of parameters 
is a special coordinate system on this space, called the spectral coordinates. Expressing isomonodromy 
transformation in those coordinates gives rise to the difference Painleve equations. 

Definition 4.1. Let L(2;) be a rational 2x2 matrix function on the Riemann sphere satisfying the following 
conditions: 

L(z)=diag{pi,p2}+-^ + -^, Li = a,bl, det L(z) = pip2 ~ f^j^^! ~ f j ■ (4.1) 

Z — Z\ Z — Z2 [Z — Zi){Z — Z2) 

In addition, put 

pifci = (Loo)ii, /02fc2 = (Loo)22, M=(Loo)2i, whcrc Lqo = - rcSoo L(z) dz = Li -|- L2 . (4.2) 

We call {pi, P2,CijC2j Zl, Z2,ki,k2, fJ.) the type of L(z). These parameters are not independent, since 

ki + k2 = trLg ^Loo = {zi - Ci) + (^2 - C2)- 

The conditions on piki correspond to fixing the formal type of the solution of the difference equation at 
infinity, and the choice of p corresponds to fixing the gauge under the global action by constant non- 
degenerate diagonal matrices. 

Definition 4.2. The spectral coordinates (7,7r) are defined by the conditions that 
• ^(7)21 = (and therefore L(z)2i = (^^^^^^i^); 

•-=8^L(7)n. 

The normalization conditions in this definition are chosen to match the formulas in [AB06] and [AB07]. 



6 



A. Dzhamay 



Factorization of Matrix Functions 



Lemma 4.1. Let L(z) be a rational 2x2 matrix function on the Riemann sphere that has the type 
{PijP2,CiX2,zi,Z2,ki,k2,lJ')- Then, in spectral coordinates, the residues of the matrix ^{z) are given by 



L2 = M 



J- Zi 
Z2 - zi 

1 - Z2 
Zl - Z2 



I {piiki + -1 - Z2) - Tl 



(7-Z2)(7^C2) 

1-zi 



(pi(fci+7-2;i)-7r(7-C2))' 



1 i 



[1 i(p2(fc2+7-^2)-^(7-Ci))], (4.3) 

I (P2(fc2 + 7 - ^1) - ^^^TFSy^)] • (4.4) 



Proof. Let L, = 



[1 6i] . Then the normaUzation at infinity results in the equations ai + a2 = 

piki = a\ax + a2a2 and ^2^:2 = ai^i + 0:2^2. The first of these equations, together with the condition 
L(7)2i = gives the expressions for aj: 



oti = 11- 



1-Zi 
Z2- Z\ 



a2 = H- 



1- Z2 
Zl - Z2 



The equation on L(7)ii is 



T / X (7-C2) , Qicii . 0:202 
L(7)ii = -K- = pi H h 



which then gives 
1 



P\{ki + 7 - ^2) - TT 



(7 - Zl) 

(7- ^2)(7- C2) 

7-2:1 



7 - 2;i 7 - -22 



(4.5) 



(4.6) 



(12 = - (pi (fci + 7 - -Zl) - 7r(7 - C2)) • (4.7) 



Using the equation L(7)iiL(7)22 — detL(7) we get L(7)22 = (7-22) ' '^'^^^ allows us to find the expres- 
sions for h: 



61 = ^ (p2(fc2 + 7 - ^2) - — (7 - Ci)) , h2 = - (p2{k2 + 7 - ^1) - — 
which completes the proof. 

Corollary 4.2. In the same gauge, the residues Mj of the inverse matrix 

Ml M2 



P1P2 (7 - zi){-/ - Ci) 



(7 - Z2) 



, (4.E 



M(2) = L(z)-i = diag{l/pi, I/P2} - 



z-Ci Z-C2 



□ 



(4.9) 



are given by 

M 7 - Ci 

Pl/'2 C2 - Cl 

A* 7 - C2 



Ml 
M2 



i(p2(fci-7 + C2) + ^(7-^i)) 



P1P2 Cl - C2 



Mp2(fci-7+Ci) + ^ 



P1P2 (7-Cl)(7-zi) 



(7-C2) 



1 ^ (pi(fc2-7 + C2) + 7r (^-gl(^y-) )] ; (4.10) 
[1 ^(pi(A;2-7 + Ci) + 7r(7-22))] . (4.11) 



Proof. Since M(2;) has the same form as L(2;), we only have to determine the type and spectral coordinates 
of yi{z) in terms of those of L(2;), 



L(z): Zl 22 Cl C2 Pi P2 ki k2 

M(z): Cl C2 Zl 22 - - -ki -k2 ^ 7 ^7~"^\^7~S^\ 

w ^1 •sz 1 z P2 ^ pi/92 ' 'r(7-22)(7-C2) 

and take the negative sign in the definitions of Mj into account 



II 7 



□ 
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CoroUciry 4.3. Using the notation f{a,b) = pi{'y — a) — Tr{j — b) we can write the formulas for L(2;) and 
M.{z) in the following symmetric form: 



Li = M 

L2 = /U 

Ml = 
Mo = 



1 - Zi 
Z2 - Zl 

I-Z2 
Zl - Z2 

M 7 - Cl 
P1P2 C2 - Cl 

M 7- C2 

P1P2 Cl - C2 



^(pifcl + ^^(^l,C2)) 



7-21 
1 

i {piki+ip{zi,C2)) 



[1 ^(p2fc2-f^(Cl,^2))], 



1 MM2-f^^v^(Ci,^2) 



(p2fcl+f^(^l,C2)) 



1 ^(Plfc2-^V(C1,^2)) 



^(Mi + f^^^(^i,C2)) 



[1 i(pifc2-^(Cl,^2))]. 



(4.12) 

(4.13) 
(4.14) 

(4.15) 



4.2 Difference Painleve V 

Consider now the isomonodromy transformation given by Ii,{z) = Bi(z): 

L(^) = Bi{z)-LoBl{z) ^ L(^) = Bl{z + l)B^(^)Lo = B'i(^)LoB^(^). 
Theorem 4.4. The type and spectral coordinates of'L{z) in terms of those of'L{z) are given by 



Mz): 
l{z): 



Zl Z2 Cl C2 Pl P2 ki h fJ- J 

zi = zi-l Z2 = Z2 Cl = Cl - 1 C2 = C2 Pl P2 ki k2 p J n 



where 



P = P 



J + 1 = Z2 + C: 



, Pl(7r - P2) 
f)2{TT - Pl) 

Pijki - Zi + C2) _^ P2{k2 - -Zl + C2 + 1) 
TT - Pl TT - P2 



^ ^ (7-g2)(7-C2) 

TTTT = plp2 — 

(7- zi)(7- Cl) 



(4.16) 
(4.17) 

(4.18) 
(4.19) 
(4.20) 



Equations (4-19-4-20) are the difference Painleve V equations of the Sakai's hierarchy [SakOlJ, first obtained 
in this setting in [AB06] (Theorem B). 

Proof. First note that 



Lcx) — G2L0 + LqGi, L(x) — (Gi + G2)Lo. 
Thus, using (3.8), (4.12), and (4.14), we get 

M = (i'<x))21 = [G5;,Lo]21 = P+{pl - P2){G:\)21 = P+{pi - P2)^ 

h[ci 



(4.21) 



M 1 + 



tt{pi - P2) 



Pli-K- P2) 



P2(7r-pi)/ P2{lT - Pl) 

Prom the uniqueness of the left and right divisors we see that B\{z + 1) = &i{z) and B2(z) = ^"62(2). 
Using the first equation, we see that GJ = Gi and so (cibi)/(biCi) = (aidi)/(diai). In particular. 



^blci = {Zl - C,)P^^^^ = dlai = {~zi - Cl) ( Pl 



7r(7- ^2)(7- C2) 
(7-^i)(7-Ci) , 



(4.22) 



which gives (4.20). Similarly, comparing the first components of the normalized vectors Ci and ai gives 
(4.19) and completes the proof. □ 
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